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To systematically describe evaporation spectra for light and heavy compound nuclei over a large 
range of excitation energies, it was necessary to consider three ingredients in the statistical model. 
Firstly, transmission coefficients or barrier penetration factors for charged-particle emission are 
typically taken from global fits to elastic-scattering data. However, such transmission coefficients 
do not reproduce the barrier region of evaporation spectra and reproduction of the data requires 
a distributions of Coulomb barriers. This is possibly associated with large fluctuations in the 
compound-nucleus shape or density profile. Secondly for heavy nuclei, an excitation-energy de- 
pendent level-density parameter is required to describe the slope of the exponential tails of these 
spectra. The level-density parameter was reduced at larger temperatures, consistent with the ex- 
pected fadeout of long-range correlation, but the strong A dependence of this effect is unexpected. 
Lastly to describe the angular-momentum dependence of the level density in light nuclei at large 
spins, the macroscopic rotational energy of the nucleus has to be reduced from the values predicted 
with the Finite-Range Liquid-Drop model. 



I. INTRODUCTION 

The statistical model of compound-nucleus (CN) decay 
is extensively used in pure and applied nuclear science. In 
many reaction scenarios one or more compound nuclei are 
formed after a nuclear collisions. Compound nuclei are 
equilibrated in their non-decay degrees of freedom and 
thus their decay is independent of how they were created. 
Statistical-model codes have been used as "afterburners" 
in many reaction-modelling programs decaying the sim- 
ulated compound nuclei produced from some initial fast 
reaction mechanism. The initial reaction could be fusion, 
spallation, fragmentation, etc. 

Accurate determination of the statistical-model pa- 
rameters at high excitation energies would give insight 
into properties of hot nucleus. The level densities are 
sensitive to the magnitude of long-range correlations as- 
sociated with collective excitations, transmission coeffi- 
cients are sensitive to the charge and mass distributions 
and fission maybe sensitive to the nuclear viscosity. Ef- 
forts to extract such information require systematic stud- 
ies of compound-nucleus decay covering a large range of 
compound-nucleus Z ,A and excitation energy. 

The modelling of spallation reactions is important 
in applications ranging from transmutation of nuclear 
waste, the design of neutron sources for condensed-mater 
studies, radiation protection around accelerators and in 
space, and for the production of rare isotopes for nu- 
clear and astrophysics experiments. The modelling of 
such reactions involves an Intra-Nuclear Cascade [l], Q or 
Quantum Molecular Dynamics code Q to simulate the 
production of the initial fast reaction products and the 
properties of the residual compound nuclei formed (Z, A, 
E* , and J joint distributions). These compound nuclei 
then de-excited with a statistical-model code which in- 
cludes evaporation and fission and possible other decay 
modes. The residual nuclei are predicted to be excited 
to large excitation energies (many hundreds of MeV) and 
therefore knowledge of the statistical-model parameters 



is needed for this energy regime. 

The final predictions of spallation modelling are sensi- 
tive to both the statistical-model parameters and those 
associated with the initial fast phase of the reaction. 
When fitting experimental spallation data, it is not al- 
ways possible to isolate the role of the statistical-model 
parameters and constrain them. Alternatively heavy-ion- 
induced complete-fusion reactions can be used to create 
compound nuclei. In complete fusion, the excitation en- 
ergy and identity of the compound nucleus are completely 
defined from conservation laws. The CN spin distribution 
can also be well constrained. The maximum spin can be 
determined from measurements of the total fusion cross 
section or, alternatively, simple one-dimensional mod- 
els are generally quite accurate above the fusion barrier. 
Thus the simple complete-fusion mechanism with no fast 
non-statistical particles and a well defined distribution of 
CN provides an opportunity to constrain the statistical- 
model parameters. 

Of course complete-fusion reactions are limited by 
preequilibrium emissions and incomplete- fusion processes 
which sets it at large bombarding energies (> 10 
MeV/ A). However large excitation energies (up to ^250 
MeV) can still be probed with complete fusion using 
more symmetric reactions. Heavy-ion-induced fusion re- 
actions, especially the more symmetric cases, emphasize 
large spins, typically larger than those probed by spalla- 
tion reactions at the same excitation energies. Therefore 
application of statistical-model parameters determined in 
fusion reactions to spallation modelling requires a good 
understanding the spin dependence of CN decay. 

The statistical model has a long history in heavy-ion 
induced fusion reactions and has been fit to a large body 
of data including fission probabilities, light-particle evap- 
oration spectra, residual Z and A distributions, gamma- 
ray multiplicities, etc. Although such data are usually fit 
within the statistical-mode framework, it has generally 
been found necessary to fine tune the statistical-model 
parameters for a particular compound nucleus or mass re- 
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gion. No statistical-model prescription exists which gives 
accurate predictions of these quantities over the entire 
table of isotopes. This work starts to address these prob- 
lems by concentrating on light-particle evaporation which 
is sensitive to the excitation energy and spin dependences 
of the nuclear level density and the transmission coeffi- 
cients for penetration of the Coulomb barriers hindering 
particle emission. 

The assumption that the decay of the compound nu- 
cleus is independent of how its was created may not al- 
ways be correct in fusion reactions. At high excitation en- 
ergies when the statistical lifetime approaches the fusion 
timescales, dynamical effects may occur which depend 
on the entrance-channel mass asymmetry. Specifically 
symmetric reaction channels are predicted to dissipate 
the entrance-channel kinetic energy more slowly and may 
start particle evaporation before the fusion dynamics is 
complete. There have been many studies of entrance- 
channel dependence of compound-nucleus decay. How- 
ever taken as a whole, no clear consistent picture has 
emerged from these studies and in a number of cases their 
conclusions are contradictory. In particular concerning 
the shapes of evaporation spectra, one should note three 
studies where a-particle spectra were measured for dif- 
ferent entrance channels, but with matched excitation- 
energy and spin distributions. Cinausero et al. found 
no entrance-channel dependence of the spectral shape 
for A —160 compound nuclei at E* —300 MeV formed 
in 86 Kr,+ 76 Ge, 16 O+ 150 Sm, and 60 Ni+ 100 Mo reactions 
01 . For £*=170 MeV 164 Yb compound nuclei formed 
in 16 0+ 148 Sm and 64 Ni+ 100 Mo reactions, Charity et al. 
noted a slight enhancement in the a-particle yield in the 
subbarrier region, otherwise the kinetic-energy spectra 
were consistent [5j. On the other hand, Liang et al. re- 
ported on entrance-channel dependences of the slope of 
the high-energy tail in £ , *=113-MeV 156 Er compound nu- 
clei formed in 12 C+ 144 Sm, 35 Cl+ 121 Sb, and 60 Ni+ 96 Mo 
reactions Q. It is difficult to reconcile there three stud- 
ies as they pertain to the same mass region. In this work 
we will ignore such effects and assume, that if they exist, 
they are small at least compared to the overall variations 
due to the mass, excitation-energy, and spin dependences 
of the statistical-model parameters 

Statistical-model parameters are extracted from com- 
parison of statistical-model calculations to experimen- 
tal data. In this work, all statistical-model calculations 
were performed with the code GEMINI++ J] writ- 
ten in the C++ language. This is a successor of the 
well known statistical-model code GEMINI [8] written 
in FORTRAN. 



tation energies and references are listed in Table Q] For 
compound nucleus with A > 150, only studies where light 
particles were detected in coincidence with evaporation 
residues were used. For the lighter systems, only inclu- 
sive spectra are available. By appropriate selection of de- 
tection angle [backward (forward) angles for normal (re- 
verse) kinematics reactions], one can isolate proton and 
a-particle spectra which are dominated by compound- 
nucleus emission though some contamination from other 
reaction processes is possible for the lowest kinetic ener- 
gies This will be discussed in more detail in Sec. fVl 

While the residue-gated spectra may be cleaner, they 
may suffer from distortions due to the limited kinematic 
acceptance of the residue detectors. For example, de- 
tection of evaporation residues at large angles enhances 
high-energy particles as these give the largest recoil kick 
to the residue enabling it to get to such angles. The 
spectra used in this work were either corrected for this 
effect in referenced studies, or, for the 160 Yb compound 
nucleus, the GEMINI-I — h simulations were gated on the 
experimental residue acceptance. 

It is important in the statistical-model calculations to 
have realistic spin distributions for the compound nuclei. 
The fusion cross section as a function of spin was assumed 
to have the form: 



(Tf us (J) = ttA 2 ^- 



(2J+1) 



(1) 



The quantity Jo can be constrained from the fusion cross 
section. This is either measured, constrained from sys- 
tematics, or obtained from the Bass model [l(| which is 
reasonably accurate for the systems under study. The 
parameter S J was varied from 2 to 10 h with increas- 
ing asymmetry of the entrance channel. However in this 
work, the sensitivity of the predicted evaporated spectra 
to this parameter is very small. 

Fission competition is also important for determin- 
ing the J values which give rise to evaporation residues. 
When available (Table I3J, fission and/or evaporation 
residue cross sections were fit by adjusting the fission 
parameter <2//a n of Sec. IIVD1 Otherwise, interpolated 
values of a//a n were used. A more detailed discussion of 
the fission parameters in GEMINI-I — h can be found in 
Ref. [TD. 



III. EVAPORATION FORMALISM 



II. DATA 

The data used in this work to constrain the statistical- 
model parameters has come from many experimen- 
tal studies covering a wide range of compound-nucleus 
masses. The compound nucleus, the reactions, the exci- 



As GEMINI++ is to be used for CN with high spins, 
the evaporation of light particles is treated with the 
Hauser-Feshbach formalism t 32| which explicitly takes 
into account the spin degrees of freedom. The partial 
decay width of a compound nucleus of excitation energy 
E* and spin Jcn for the evaporation of particle i is 
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TABLE I. Experimental data used in this work indicating the compound nucleus (CN), the beam energy Ebeam, the excitation 
energy E* , the fusion reaction, the evaporation spectra measured (n,p,a), the values of Jo defining the angular-momentum 
distribution of Eq. ([TJ) . The first listed reference refer to the study that measured the kinetic energy spectra. The a references 
refers to measurements of the fission and residues cross sections used to constrain Jo and the fission probability. 



CN 


^beam 

[MeV] 


E* 
[MeV] 


reaction 


ref. 


spectra 


a refs. 


Jo 


b9 Cu 


100 


58 


3a S+ ay AI 


[12] 


a 


[13-16] 


27^, 




105 


60 




[17] 


a 


[13 16] 


30f 




130 


72 




[12] 


a 


[13-ro] 


34^ 




140 


77 




[12] 


a 


[13 16] 


38^ 




150 


82 




[12] 


a 


[13-i6] 


39^ 


67 Ga 


214 


110 




[17] 


a 


[13-ro] 


45^ 


187 


90 


40 Ar+ 27 Al 


[18] 


p,a 




46^ 




670 


127 


55 Mn+ 12 C 


[19] 


p,a 




42^ 




280 


127 


40 Ar+ 27 Al 


[19] 


p,a 




54^ 


96 Ru 


180 


113 


32 S+ 64 Ni 


[20] 


p,a 




69^ 


106 Cd 


160 


99 


32 S+ 74 Ge 


[21] 


p,a 




683 




99 


291 




[21] 


p,a 




83^ 




99 


291 




[21] 


p,a 




89^ 




99 


291 




[21] 


p,a 




893 


117 Te 


81 


71 


14 N+ 103 Rh 


[22] 


p,a 




403 




146 


71 


40 Ar+ 77 Se 


[23] 


p,a 




523 




121 


106 


14 N+ 103 Rh 


[22] 


p,a 




53 3 


156 Er 


142 


113 


12 C+ 144 Sm 


m 


p,a 


[6] 


543 




218 


113 


35 Cl+ 121 Sb 


[6] 


P 




>863 




333 


113 


60 Ni+ 96 Zr 


_ 


p,a 


[24] 


>903 


160 yb 


300 


91 


60 Ni+ 100 Mo 


[25] 


n,p,a 


[25] 


>903 




360 


129 




[25] 


n,p,a 


[25] 


>903 




420 


166 




[25] 


n,p,a. 


[25] 


>903 




480 


204 




[25] 


n,p,a 


[25] 


>9o3 




546 


245 




[25] 


n,p,a 


[25] 


> 9(7 


193rpj 


145 


65 


28 Si+ 160 Ho 


[26] 


p,a 


[26] 


46^ 




166 


83 




[26] 


p,a 


[26] 


683 




193 


106 




[26] 


p,a 


[26] 


803 




216 


125 




[26] 


p,a 


[26] 


963 


200 Pb 


121 


86 


i9 F+ isi Ta 


[27] 


p,a 


[27] 


563 




154 


116 




[27] 


p,a 


[27] 


72 T 




179 


139 




[27] 


p,a 


[27] 


833 
89^ 




195 


153 




[27] 


p,a 


[27] 


224 Th 


114 


59 


16 0+ 208p b 


[26] 


p,a 


[28-30] 


5l3 


224 Th 


138 


82 




[26] 


p,a 


[28-30J 


673 



a The Jo is constrained from the measured fusion cross section. 

b The Jo values were obtained from the referenced studied where the fusion cross section was estimated from the systematics of Ref. |3l| . 
c The Jo values were estimated from the Bass model [ld |. 

d The Jo vlaues are large and the residue cross section is determined solely by fission competition. Fission parameters were adjusted to 
reproduce measured evaporation residues. 



barrier penetration factor, and pd and pcn are the level 
densities of the daughter and compound nucleus, respec- 
tively. The summations include all angular momentum 
couplings between the initial and final states. In GEM- 
INIH — h, the Hauser-Feshbach formalism is implemented 
for the n, p, d , t, 3 He, a, 6 He, 6 " 8 Li, and 7 - 10 Be 
channels. However in this work, we will just compare 
predicted kinetic-energy spectra to experimental results 
for the p, a, and occasional n channels. GEMINIH — h 
also allows for intermediate-mass fragment emission fol- 
low the formalism of Moretto [33|. However, these de- 
cay modes are not very important for calculations of this 



1 r 00 

T % {E*,J CN ) = — / ds 2 

2<KpcN {E*, Jcn) J f^ Q 

JCN + Jd J+Si 

53 53 T t (e)p d (E*-B i -s,J d ) (2) 

J=\JCN-Jd\ l=\JSi\ 

where Jd is the spin of the daughter nucleus, Si, J, and 
£, are the spin, total and orbital angular momenta of 
the evaporated particle, e and J?.; are is its kinetic and 
separation energies, Tg is its transmission coefficient or 
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work. 

The nuclear level density is often approximated by the 
Fermi-gas form [34j derived for a spherical nucleus in the 
independent-particle model with constant single-particle 
level densities; 



Pfg(E*,J) 



(2 J + 1) 



24V2a 1 /4 C/5/4 CT 3 
S = 2VaU 



exp(S) 



(3) 
(4) 



where S is the nuclear entropy and the level-density pa- 
rameter is 



(5) 



Here g n (£p) and g p {s p F ) are the neutron and proton 
single-particle level densities at their respective Fermi en- 
ergies and 

U = E* — E rot (J), E rot = J{J J 1)H2 , (6) 



(7) 



The quantity I r i g is the moment of inertia of a rigid body 
with the same density distribution as the nucleus and T 
is the nuclear temperature; 



1 

T 



dS_ 

dU 



(8) 



The quantity U can be interpreted as a thermal excita- 
tion, after the rotational energy of the nucleus is removed. 

At large angular momenta, macroscopic models of 
the nucleus such as the Rotating Liquid-Drop Model 
(RLDM) [35| and Sierk's Yukawa-plus-exponential finite- 
range calculations [IB] predict the nuclear shape distorts 
to accommodate the centrifugal forces. Many implemen- 
tations of the statistical model including GEMINI-) — h, 
generalize Eq.0 by the replacing E rot (J), the rotational 
energy of a spherical nucleus of fixed moment of inertia, 
with E yarst (J), the deformation-plus- rotational energy 
predicted by these macroscopic models where the defor- 
mation increased with spin. In GEMINI -| — h, the Sierk 
predictions of E yrast {J) are used for all but the lightest 
compound nuclei (see Sec. W\. 

The shape of the kinetic-energy spectra of evaporated 
particle is thus sensitive to three ingredients. 

• The magnitude of the level-density parameter and 
its excitation-energy dependence 

• The transmission coefficients Ti (e) . 

• The angular- momentum dependence of E yrast (J). 

The level-density parameter defines the slope of the 
exponential tail of the evaporation spectrum while the 
transmission coefficients define the shape in the Coulomb 
barrier region and the effects of these two ingredients are 



easily isolated when comparing to data. The angular- 
momentum dependence of E yrast {J) is most important 
in light nuclei where the moments of inertia are small 
and thus E yrast (J) rises rapidly with spin. In particu- 
lar E yrast has a strong influence on the heavier fragments 
such as a particles which can remove large amounts spins. 
For these particles, the functional form of E yrast (J) can 
make significant modifications to the predicted shape of 
the evaporation spectrum in the exponential tail and even 
in the Coulomb barrier region. The effect of E yrast can 
be disentangled from the effects of the level-density pa- 
rameter and the transmission coefficients by comparing 
data for a lighter particle such as a proton to that for a 
heavier particle such as an a particle. 

In the following three sections, the parametrization of 
these three ingredients needed to describe experiment 
data will be described. We will discuss light and heavy 
systems separately. 



IV. HEAVY COMPOUND NUCLEI 

Let us start by concentrating on the heavier compound 
nuclei with A >I50 for which the evaporation spectra are 
shown in Figs.[T]to[5j These data sets were all obtained 
with a coincidence requirement of a detected evapora- 
tion residue. We will first consider which transmission 
coefficients and level-densities allow us to reproduced the 
shape of the experimental spectra. Predicted spectra in 
these figures will be normalized to give the same peak 
differential multiplicity in order to concentrate of the re- 
production the spectral shapes. Subsequently we will re- 
turn to consider how well one can reproduce the absolute 
multiplicities of evaporated protons and a particles. 



A. Transmission Coefficients 

The evaporation formalism is justified on the condition 
of detailed balance. The evaporation rate of an isolated 
compound nucleus is assumed identical to the emission 
rate of such a nucleus in equilibrium with a gas of the 
evaporated particles. In equilibrium there is a balance 
between the emission and the inverse, absorption rates 
of that particle and thus the transmission coefficients or 
barrier penetration probabilities should be identical to 
those for the inverse absorption process. 

Transmission coefficients have traditionally been ob- 
tained from the inverse reaction using optical-model 
parameters obtained from global optical-model fits to 
elastic-scattering data. There are t wo p roblems with this 
approach. First, Alexander et al. [3 71 ] have pointed out 
that such transmission coefficients contains the effects 
of transparency in the inverse reaction which is not ap- 
propriate in evaporation. Instead is was suggested that 
the real optical-model potentials should still be used, but 
to ensure full absorption, the incoming-wave boundary- 
condition (IWBC) model [38| be used to calculate T(. In 
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V:<a) E*=153 MeV 
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FIG. 1. (Color online) Center-of-mass kinetic-energy spectra 
of a particles and protons detected in coincidence with evap- 
oration residues formed in 16 O+ 208 Pb reactions. Experimen- 
tal results (data points) are shown for the indicated excitation 
energies of the 224 Th compound nuclei. The curves show spec- 
tra predicted with GEMINlH — (- code and normalized to the 
same peak height as the experimental data. The solid curves 
(the default calculations of the code) were obtained with the 
excitation-dependent level-density parameter and with distri- 
butions of Coulomb barriers. The short-dashed curves in- 
dicated the results obtained using a single Coulomb barrier 
and the long-dashed curves are associated with an excitation- 
independent 2=A/7.3 MeV -1 level-density parameter. 



GEMINI++, global optical- model potentials were ob- 
tained from Refs. [39l - l45| . The difference between IWBC 
and optical-model transmission coefficients is only im- 
portant for neutrons, protons, deuterons, tritons, and 
3 He particles as other particles experience strong absorp- 
tion inside the Coulomb barrier. Due to transparency, 
optical-model transmission coefficients for nucleons do 
not approach unity for energies well above the barrier 
as is the case the for IWBC values. However, the dif- 
ference between IWBC and standard optical-model Tt 
values is not that large and it is difficult to differentiate 
them based on experimental data due to uncertainties 
in other statistical-model parameters. Comparisons of 
statistical-model predictions with IBWC and standard 
optical-model value of Tt are made in Ref. [2(J [4(| where 
the biggest differences are associated with deuteron and 
triton spectra; 

The more important problem with the traditional 
transmission coefficients is that they are not associated 
with the inverse reaction. The true inverse process to 
evaporation is the absorption of the particle by a hot, 
rotating target nuclei which is impossible to measure 
experimentally. This is highlighted by the fact that 
IWBC and optical-model transmission coefficients fail to 
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FIG. 2. (Color online) As in Fig. Q] but now for 200 Pb com- 
pound nuclei formed in 19 F+ 181 Ta reactions. 
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FIG. 3. (Color online) As in Fig. Q] but now for 193 T1 com- 
pound nuclei formed in 32 Si+ 160 Ho reactions. 
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20 
e [MeV] 

FIG. 4. (Color online) As in Fig. [Q but now for 160 Yb com- 
pound nuclei formed in 60 Ni+ 100 Mo reactions with neutron 
spectra also included. 
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FIG. 5. (Color online) As in Fig. [T] but now for 156 Er com- 
pound nuclei formed in the three indicated reactions. 



reproduce the shape of the low-energy or "sub-barrier" 
region of the spectra of a and other heavier particles 
[1 S, M, HE Wi EzHH. We illustrate this in Figs. Q] 
to [5] where statistical-model predictions obtained with 
GEMINIH — h using the IWBC transmission coefficients, 
indicated by the short-dashed curves, are compared to 
experimental data. The level-density prescription used 
in these calculations will be described in the following 
sections and Sierk's values of E yrast were used. For a 
particles emitted from these heavier systems, the rela- 
tive yield in the "sub-barrier" region is clearly underpre- 
dicted. 

Some studies have attempted to reproduce such data 
by reducing the Coulomb barrier, for example by allow- 
ing an extended radial-profile of a spherical nucleus [5(| ■ 
However, a simple reduction in the barrier, just shifts 
the kinetic-energy spectrum down in energy. The exper- 
imental a-particle spectra have more rounded maxima 
than predictions with such barriers. This is illustrated in 
Fig. [5] where the a spectrum measured for E*—120 MeV 
193 T1 compound nuclei formed in 28 Si+ 160 Ho reactions is 
compared to a number of calculations. The solid curve 
is again the prediction with the standard IWBC trans- 
mission coefficients. For the short-dashed curve, the 
Coulomb barrier was decreased by increasing the radius 
parameter of the nuclear potential by Sr from its orig- 
inal value of Rq in the global optical-model potential. 
The value of Sr is temperature dependent and is given 
later. With the reduced barrier, there is a predicted in- 
creased in the yield at lower energies but the yield starts 
dropping too early with energy and doesn't reproduce 
the width of the experimental distribution. For interest 
sake, the spectrum predicted with = Rq — Sr is indicated 
by the long-dashed curve. Although decreasing the level- 
density parameter will increase the predicted width of 
the spectrum, the exponential slope of the experimen- 
tal spectrum is already reproduced for E c . m . >27 MeV 
by all the curves. It is clear that if one considered a 
distribution of radius parameters, one could increase the 
predicted width of the a-particle spectrum. This conclu- 
sion was also found for evaporated Li and Be particles 

A distribution could arise from a static nuclear defor- 
mation if evaporation is averaged over the nuclear sur- 
face [Hj]. Alternatively, the origin of this distribution 
may have contributions from compound-nucleus thermal 
shape fluctuations [H, H3] and/or fluctuation in the dif- 
fuseness of the nuclear surface or nuclear size. 

If the fluctuations are thermally induced then we ex- 
pect, to first order, their variance to be proportional to 
temperature. In GEMINI-I — h, a simple scheme was im- 
plemented to incorporate the effects of barrier distribu- 
tions. The transmission coefficients were calculated as 



Tt{e) 



T*°- Sr (e) + lf° (e) + T e Ro+Sr (e) 



(9) 



which is the average of three IWBC transmission coeffi- 
cients calculated with three different radius parameters 
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FIG. 6. (Color online) Comparison of the experimental a- 
particle evaporation spectrum (data points) measured in the 
2S Si+ 160 Ho reaction producing 193 T1 compound nuclei at 
E* = 126 MeV to GEMINI++ predictions. The solid curve 
was obtained with standard IWBC transmission coefficients, 
while the short- and long-dashed curves were obtained by in- 
creasing and decreasing the radius parameter of the nuclear 
potential, respectively, (see text). The curves have been nor- 
malized to the same peak height as the experimental data. 

of the nuclear potential. It was assumed 

8r = wVT (10) 

consistent with thermal fluctuations where the value of 
the parameter w =1.0 fm was obtained from fits to ex- 
periment data and T is the nuclear temperature of the 
daughter nucleus as defined in Eq. ((8]). An example 
of these transmission coefficients is shown in Fig. [7] for 
a+ 193 Tl with I = at T=3 MeV. The dashed curves 
show three transmission coefficients associated with the 
three radii in Eq. (J5J) and the solid curve is the final re- 
sult, the average of the three dashed curves. The more 
gradual rise of the transmission with kinetic energy gives 
rise to a broader peak in the predicted a-particle spectra. 

Results obtained with this prescription are indicated 
by the solid curves in Figs. [1] to [5] and generally reproduce 
the a particle data quite well. 

Because of their lower absolute Coulomb barriers, the 
effect of the distribution is much less for protons and is 
practically absent for neutrons. However, the agreement 
for protons is generally improved. 

One should note that the magnitudes of the fluctua- 
tions are very large. For a temperature of T=3 MeV, 6r 
is ~ 25% of the nuclear radius for A=160. For ellipsoidal 
shape fluctuations in Ref. [5lj , the full width at half max- 
imum of the Coulomb barrier distributions was predicted 
to be only ~7%. This suggests that either higher-order 
shape fluctuations are required or the fluctuations are 
associated with density profile. 

The effects of the barrier distributions is to increase 
the width of the kinetic-energy window around the bar- 
rier where the transmission coefficients change signifi- 
cantly. For example in Fig. UJ the transmission coefficient 



FIG. 7. (Color online) Transmission coefficents for q+ 193 T1 
at 1=0. The dashed curves shows the three transmission co- 
efficents which different nuclear radii which are averaged in 
Eq. (J2J) and the sold curve is the result. 



changed from 10% to 90% over an interval of 4.5 MeV for 
IWBC calculation [Tf°(e)]. However, with Eq. ©, this 
increased to 9.2 MeV. An alternative way of increasing 
the width of this window would be to make the radial 
width of the barrier narrower. Narrow barriers allow for 
more tunneling and enhance the transmission just below 
the barrier and also decrease it just above the barrier. 
However it is difficult to see how the barrier could be 
made significantly narrower as the decrease in the poten- 
tial at large distances is dictated by the Coulomb poten- 
tial which falls off slowly. Thus barrier distributions are 
the most likely explanation. 



B. Level-Density Parameter 

The slope of the exponential tail of the kinetic-energy 
spectrum gives sensitivity to the nuclear temperature T 
[Eq. ©]. The temperature is dependent on the rate of 
change of the level density, but not its absolute value. 

The Fermi-gas level density prescription of Sec. lIIII can 
be further refined by including the pairing interaction 
[55l l56| . For the spin and excitation-energy region of 
interest in this work, the pairing gap has vanished and we 
can use a back-shifted Fermi-Gas formula by substituting 
the following definition of the thermal excitation energy 

U = E*-E yrast (J)+6P (11) 

where SP is the pairing correction to the empirical mass 
formula. 

At low excitation energies, the absolute level den- 
sity can be measured via neutron-resonance counting. 
The level-density parameters extracted from such data 
in Ref. [57| . using the back-shifted Fermi-gas formula, 
are plotted in Fig. [5] The level-density parameter has 
strong fluctuations due to shell effects which can be 
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FIG. 8. (Color online) Mass dependence of level-density 
parameters. Experimental points from neutron-resonance 
counting are shown as the filled, square data points. The 
open circles are fits obtained using Eq. (|12p . 



parametrized as [58| 
a(U) = a 



SW 

l-h(U/r, + J/J v ) — 



(12) 



where SW is the shell correction to the liquid-drop mass 
and a is a smoothed level-density parameter. With 
h(x) — tanh(x) we obtain a best fit (open-circular 
points) to the experimental data with r\= 19 MeV and 
a = A/7.3 MeV" 1 . 

The angular-momentum dependence of h(x) is irrele- 
vant for neutron resonances which are S-wave in nature. 
However for fusion reactions, it was decided to include 
a fading out of shell effects which spin. Although at 
high spins and low values of U, shell corrections are still 
important, the configuration of the nucleus has changed 
from the ground state and the use of the ground-state 
shell correction is wrong. Rather than use an incorrect 
shell correction, it was decided to use no correction at 
all. The parameter was set to 50 h. 

The above prescription for the fadeout of shell and 
pairing corrections is used in all GEMINI-I — h calcula- 
tions with separation energies Bi, nuclear masses, shell 
SW and pairing SP corrections obtained from the tabu- 
lations of Mollcr et al. . 

Predicted kinetic-energy spectra obtained using these 
pairing and shell modified Fermi-gas level density pre- 
scription are shown as the long-dashed curves in Figs. Q] 
to [5j They significantly underestimate the yield in the 
exponential tails for the heavier systems. This disagree- 
ment gets worse with both increasing compound-nucleus 
mass and excitation energy. These results suggest that a 
excitation-dependent value of a is needed. 

The value of the smoothed level-density parameter a 
used in these calculations is large compared to estimates 
from the independent-particle model of a — A/ 10- A/ 11 
MeV -1 [60, 61] and the difference has been attributed to 



correlations. In particular, it is the long-range correla- 
tions associated with coupling of nucleon single-particle 
degrees of freedom to low-lying collective modes and gi- 
ant resonances which are most important. 

It has been proposed that long-range correlations mod- 
ify the Fermi-gas level density in two wa ys. The first of 
these is called collective enhancement [62ll63l] . For exam- 
ple if we have a deformed nucleus, then for each single- 
particle configuration, one can consider collective rota- 
tions. In additions, both spherical and deformed nuclei 
can have collective vibrational motions. These collective 
motions give rise to rotational and vibrational bands en- 
hancing the level density above the single-particle value, 



p{E*) = K coll (E*)p FG (E*) 



(13) 



where K co u is the collective enhancement factor. 

Long-range correlations, and to a lesser extent also 
short-range correlations, cause an enhancement of the 
sing le-particle level densities g n {e F ) and g p (e p F ) in Eq. (O 
64] which leads to an enhancement in a. This enhance- 
ment is counterbalanced by the effect of nonlocality. In 
fact without the correlations, we would expect smaller 
level-density parameters than the predicted a = A/10- 
A/ll MeV -1 values due to the unbalanced effect of non- 
locality. As U increases, long-range correlations are ex- 
pected to wash out giving rise to both a disappearance 
of collective enhancement [K co u — > 1) and a reduction in 
the level-density parameter itself [62, 63, 65]. 

In this work, we interpret level densities through the 
Fermi-gas formula, i.e., take Eq. Q as correct by defini- 
tion, but use an effective level-density parameter a e // 
that is enhanced above the single-particle estimate of 
Eq. ([5]) and decreases with excitation energy due to the 
fade out of these long-range correlations, i.e., 

p(E*) = p FG (E*,a eff ) = K coll (E*)p FG (E*,a). (14) 

At low energies, a e // is set to the value of A/7.3 MeV -1 
to be consistent with the counting of neutron resonances. 

We have parametrized its excitation-energy depen- 
dence by 



a-eff (U) 



A 



koo (koo fco)cxpf fcjlfc A 



(15) 



where ko=7.3 MeV and the asymptotic value at high ex- 
citation energy is a e ff=A/koo- The parameter k defines 
how fast the long-range correlations wash out with ex- 
citation energy. This expression is expected to be valid 
only to moderately high excitation energies where expan- 
sion and increases in the surface diffuseness [65l . |66| are 
not significant. 

Experimental evidence for an excitation-energy depen- 
dence of a e ff was found in the A ~160 region; mea- 
surements of light-particle evaporation spectra (ra, p, a) 
with excitation energies ranging from 50 to 250 MeV 
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9, 25] show clear evidence of a departure from a con- 
stant value of a e f f with the data being reproduced by 
the parametrization 



eff (U) 



A 



k + kU/A 



(16) 



when ko=7 MeV and k—1.3 McV. This equation is just 
a lower-order approximation of Eq. (1151) . From an ex- 
amination of other studies on evaporation spectra, it is 
apparent that there is a strong A dependence of k. Neb- 
bia et al. [2l| find no deviation from a constant a e // 
value for the 106 Cd CN with excitation energies up to 
291 MeV [2l|. Whereas for heavier systems, larger val- 
ues of k are deduced; values of k=2-3 were found for 
A -200 (E* <150 MeV) [H, H3 and k=8.5 for ,4=224 
(E* <90) [H with fc =8 MeV. 

In this work, we have made a systematic study of the A 
dependence of k by fitting the evaporation spectra with 
Eq. (fT5j) . At the excitation energies studied, we cannot 
constrain the value of koo and it was set to 12 MeV. The 
fitted values of k obtained from reproducing the evapora- 
tion spectra in Figs. [T] to 0] are plotted verses A in Fig. |H1 
For a single compound nucleus, the values of k obtained 
from fitting the proton and a-spectra were similar though 
not always identical and the error bars in Fig. [9] reflect 
this range of n values. 

In addition to these data points, Fig. [3] gives some lim- 
its for k obtained from 117 Te and 106 Cd compound nuclei. 
These data are in fact consistent with k=0 and will be 
discussed in more detail in Sec. IVl 

Figure[S] is a log plot and it indicates that k increases 
very rapidly with mass number. Although we do not 
have enough data points to determine this dependence 
in detail, we have fitted it with the exponential function 
shown by the solid line in this figure and given by 



k(A) = 0.00517 exp(0.0345^4). 



(17) 



The excitation-energy dependence of the level-density 
parameter associated with this dependence is illustrated 
in Fig. [TU] for the indicated A values. The excitation 
dependence is very strong for the heaviest compound nu- 
clei, but below A <100, there is very little dependence. 
Statistical-model calculations preformed with this depen- 
dence are indicated by the solid curves in Figs. Q] to [5] 
They reproduce the data much better than a constant 
a e ff=A/7.3 MeV" 1 , though they are not perfect. 

The individual fits to each reactions (not shown) are 
slightly better but quite similar. In the similar mass 193 T1 
and 200 Pb systems, the tails of the a-particles spectra are 
under and over predicted, respectively. This could just 
be an artifact due to small experimental errors in the two 
studies or may reflect an asymmetry (N — Z)/A depen- 
dence of k or even an entrance channel effect. Also the 
224 Th data clearly suffer from large statistical errors due 
to the very small residue cross sections. Further system- 
atic measurements of a large number of compound nuclei 
with the same experimental apparatus would help resolve 
these issues. 




FIG. 9. (Color online) Values of n in Eq. (|15[) obtained from 
fitting evaporation spectra. The solid line shows a smooth 
approximation used to calculate evaporation spectra and ER 
excitation functions. The dashed curve shows k values ex- 
tracted from the predictions of Ref . [65[ . 
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FIG. 10. (Color online) Excitation-energy dependence of the 
smoothed level-density parameter obtained in this work for 
the indicated A values. 



More sophisticated calculations of nuclear level density 
have been obtained within the Shell-Model Monte Carlo 
method but only for light nuclei such as 56 Fe have cal- 
culations been extended to high excitation energies [67j . 
These calculated level densities can be fit with a constant 
level-density parameter of value A/9. 5 MeV -1 . This ba- 
sically consistent with the results of this work in that the 
level density of light nuclei has a Fermi-gas form (a e ff 
independent of U), however the value of A/9.5 MeV -1 is 
a little smaller than the value A/7.3 MeV -1 used in this 
work. 



C. Multiplicities and Cross Sections 

So far we have only considered the shapes of the kinetic 
energy spectra. It is also important to determined the ac- 
curacy to which the absolute yields of evaporated parti- 
cles can be predicted. For the 156 Er, 160 Yb, 193 T1, 200 Pb, 
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FIG. 11. (Color online) Comparison of experimental and pre- 
dicted proton multiplicities from the indicated compound nu- 
clei. To aid in viewing, the data have been scaled by the 
indicated factors. The solid curves were obtained with the 
excitation-dependent level-density parameter and with dis- 
tribution of Coulomb barriers. The dashed curves shows 
the prediction with single Coulomb barriers and a constant 
a eff = ,4/7.3 MeV" 1 . 



and 224 Th compound nuclei for which light particles were 
detected in coincidence with evaporation residues, the 
predicted multiplicities are compared to the experimen- 
tal proton and a-particles values in Figs. QT] and [121 To 
separate the data from the different systems, the multi- 
plicities were scaled by the indicated amounts. The solid 
curves in both figures show calculations with the default 
setting of the code, i.e., distribution of Coulomb barri- 
ers and an excitation-dependent level-density parameter 
a e ff. They reproduce the a-particle data quite well. For 
protons, the 160 Yb and 193 T1 data are well reproduced, 
while the other systems underpredict the multiplicities 
by up to a factor of 2. 

It is difficult to understand how a better overall repro- 
duction of the experimental proton multiplicities can be 
obtained for A —160. For example the 156 Er and 160 Yb 
compound nuclei have similar Z and A values, are both 
produced in Ni induced reactions and thus explore sim- 
ilar spin distributions. The protons are predicted to be 
emitted at large excitation energies where shell and pair- 
ing effects are expected to be washed out. Modifications 
to GEMINI-I — h that increase the proton multiplicity for 
the 156 Er system will also increase the multiplicities for 
the 160 Yb system in disagreement with the experimen- 
tal data. One should consider whether the inability to 
simultaneously fits these two systems is an experimental 
problem. 

These predicted multiplicities are quite sensitive to the 
level-density and Coulomb barrier prescription. To illus- 



FIG. 12. (Color online) As for Fig. [TT] but now for a-particle 
multiplicities. 



trate this, the dashed curves in Fig. [Til and [T2l were ob- 
tained with a constant a e // = A/7. 3 MeV -1 and with 
the IWBC transmission coefficients for a single Coulomb 
barrier. For a particles, this results in a large decrease of 
the multiplicities by a factor of 3 to 10. Clearly the level 
density and Coulomb barrier distribution are important 
to correctly predict these multiplicities. For protons we 
are somewhat less sensitive to these ingredients. 



D. Consequence for Fission 

Although this work is not focused on the fission prob- 
ability, it is interesting to determine the consequences 
of the parametrizations in the preceding sections on the 
fission probability. Fission was first incorporated into 
the statistical model by Bohr and Wheeler using the 
transition-state formalism first introduced to calculate 
chemical reaction rates. The Bohr- Wheeler decay width 
[H is 



T BW (E*,J) = — — / p s (E* - B f (J) - e, J)de 

PCN\E ,J) J 

(18) 

where Bf(J) is the spin-dependent fission barrier, and ps 
is the level density at the transition state, i.e., the saddle- 
point configuration. The variable e is the kinetic energy 
in the fission degree of freedom at the saddle-point. Later 
in a one-dimensional diffusion model, Kramers [69| de- 
rived a formula similar to this with a different factor be- 
fore the integral. For large viscosity, the decay width 
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fk 



,{e*,j) = f k vf w (E*,j), 



(J_\ _ 2 



(19) 
(20) 



where 7 is the magnitude of the viscosity and u is the cur- 
vature of the potential energy at the saddle-point. The 
Kramers factor //. scaling the Bohr- Wheeler width is less 
than unity and is hard to extract experimentally due to 
the much larger uncertainty associated with the fission 
barrier and the level-density parameter. 

The fission decay width has also been suggested to be 
transient [70j, i.e. initially zero and then rising to the 
quasi-stationary value of Kramers. This idea has helped 
to explain the larger number of neutrons emitted before 
the scission point is attained |7lj . During the transient 
time which can also be thought as a fission delay, any 
light- particle evaporation will lower the excitation energy 
and spin of the decaying nucleus and subsequently may 
reduce its fission probability. 

However, there is some controversy as to whether tran- 
sient fission decay widths are needed to explain exper- 
imental fission probabilities. A number of theoretical 
studies reproduce experimental fission probabilities and 
pre-scission neutron multiplicities with transient fission 
widths [72|,[73j]. The viscosity which determined the tran- 
sient time scale was found to increase with the mass in 
these studies. Transient fission has also been invoked 
to explain the unexpectedly large number of evaporation 
residues measured in the very fissile 216 Th compound sys- 
tem formed in 32 S+ 184 W reactions [74]. Alternatively 
other studies have reproduced fission probabilities [75[ 
and both prescission neutron multiplicities and fission 
probabilities with no transient effects [76|. Similarly, in 
very-high-excitation-energy data obtained with 2.5-GeV 
proton induced spallation reactions, no transients were 
needed in reproducing the measured fission yields fn\ . 

In this work, we will not try and answer all these un- 
certainties pertaining to fission, but will investigate how 
the excitation-dependent level-density parameter affects 
the fission probability. The fission decay width will be 
taken from the Bohr- Wheeler formalism. Let us assume 
that the level-density parameter for the saddle-point and 
ground-state configurations are identical apart from a 
scaling factor a/ /a n which accounts for the increased sur- 
face area of the former [7q |. Fission decay widths were 
calculated using the angular-momentum-dependent fis- 
sion barriers of Sierk [36|. For 200 Pb, 216 Th, 224 Th, and 
224 Ra compound nuclei formed in the reactions listed in 
Table HH both ER and fission excitation functions have 
been measured allowing us to determined the fusion cross 
section and thus constrain the CN spin distributions. 
Evaporation-residue excitation functions were calculated 
with the exponential dependence of k in Fig. [9] and some 
final adjustment was made with the parameter a//a n in 
order to reproduce the experimental data. The results, 
shown by the solid curves in Fig. 1131 reproduce the data 
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FIG. 13. (Color online) Evaporation-residue excitation func- 
tions for the indicated compound nuclei. The data points are 
published experimental results and the short-dashed, long- 
dashed, and solid curves were calculated with a e //= A/7.3, 
a=A/\\ MeV -1 and Eq. ([15]), respectively. 



TABLE II. Experimental data used in Fig 1131 are listed with 
the compound nucleus, reaction, references and a / /a n values 
used in the GEMINI+ + calculations. 



CN 


reaction 


ref. 


a f /a n 


2UU Pb 


i« F+ m Ta 


[27, 79, 80] 


1.04 


216 Ra 


19 F+ 197 Au 


W 


1.04 


2i6 Th 


32g + 184 w 


[74, 82] 


1.06 


224 Th 


16 + 208p b 


[26, 28-30] 


1.035 



quite well and the fitted a//a„ values, which are all sim- 
ilar in magnitude, are listed in Table [TTT For comparison, 
the short-dashed curves show the results obtained with 
a constant a e //=A/7.3 MeV -1 . The U dependence of 
a e ff gives rise to an enhancement of the predicted ER 
yield which is most pronounced for the heavier systems 
and the higher excitation energies. However, for the en- 
ergy regime where there is significant enhancement, the 
fission cross sections are orders of magnitude larger and 
even with this enhancement, ER survival is still a rare 
process. 

The calculations with the excitation-dependent values 
of a e ff have higher nuclear temperatures than the A/7.3 
MeV -1 calculation. Larger temperatures enhance rare 
decay modes and these rare decay modes are the evapo- 
ration channels in these very fissile nuclei. This is illus- 
trated by the long-dashed curves which are calculations 
with a constant a e f f—A/11 MeV -1 where the tempera- 
tures are 20% larger than for a e ff— A/7.3 MeV -1 . These 
curves also show enhanced evaporation residue yields, 
but the excitation-energy dependence is not as well de- 
scribed as by the solid curves with the excitation-energy 
dependence. For the 216 Th system of Fig. [TBT ah Back 
et al. using calculations with constant a e // concluded 
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that the statistical model was not able to reproduced 
the data and thus deduced that there must be fission 
transients 30]. However, it is now clear that with an 
excitation-energy dependent a e ff, this conclusion is no 
longer valid. This suggests possibility a reduced role for 
fission transients in determining the fission probability. 

It should be noted that the ability of these calcula- 
tions to reproduce the evaporation-residue cross sections 
depends on the assumed excitation-energy dependence 
of a e ff. For the lighter 160 Yb system of Fig. HI the 
excitation-energy dependence is rather well established 
[25j . A larger range of compound- nucleus excitation ener- 
gies were probed (Table HJ and neutron evaporation spec- 
tra were also measured. Charged particles are typically 
emitted early the decay chain and probe higher excitation 
energies whereas neutrons are emitted at all decay stages 
and give information more on the average temperature. 
Reproduction of both charged-particle and neutron spec- 
tra required an excitation-energy dependence of a e / / for 
160 Yb. Subsequently these a e // values were found con- 
sistent with data from the similar-mass 178 Hf compound 
nuclei at even lower excitation energies Q. By contrast 
only charge-particle spectra were measured for the 224 Th 
compound nucleus in Ref. (2(| and at just two excitation 
energies separated by ^20 MeV. It also was possible to fit 
these spectra with a constant a e ff=A/15 MeV [26|. Al- 
though a constant value is unlikely given the larger values 
derived from counting neutron resonances, it is clear that 
for this heavy nucleus, the excitation-dependence of a e // 
is not well constrained from the present experimental 
data. Clearly further experimental studies of this point 
would be useful in understanding the fission in these very 
heavier systems. Also it should be noted that for A ~220, 
quasifission also competes with fusion reactions at the 
lower I waves associated with evaporation-residue pro- 
duction for entrance channels with 19 F projectiles and 
heavier [8ll . l83j . This suggests that somewhat smaller 
values of aj/a n are associated with the 19 F+ 197 Au and 
32 S+ 184 W reactions than those of Table HU 

Finally it is of interest to consider the relevance of this 
work to the production of superheavy elements. Of par- 
ticular interest are "hot" fusion reactions which have pro- 
duced the heaviest elements to date [13, HE] ■ Based on 
an extrapolation of K to the A =277-294 region we would 
expect significantly enhanced temperatures for the CN 
excitation energies of ^35 MeV produced in these reac- 
tions. Therefore, this effect may also contribute to an 
enhanced yield of superheavy elements in these hot fu- 
sion reactions. Clearly, more studies are also needed in 
this area. 



E. Thermal Properties of Nuclei 

The thermal properties of nuclei can be inferred from 
the level density [34|. Figures [TWa) and [T4T b) show the 
excitation-energy dependence of S and T plotted in a 
manner that the mass dependence would disappear for 
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FIG. 14. (Color online) Excitation-energy dependence of (a) 
the nuclear temperature and (b) the entropy deduced in this 
work. 



an energy- independent a e // oc A. The curves for dif- 
ferent masses are only plotted up to the maximum U 
sampled in the experiments. We see a small mass depen- 
dence of S/A, but a larger dependence for the tempera- 
ture. For a given U/A, we see smaller values of S/A and 
larger temperatures for the heavier systems. The larger 
temperatures are responsible for the stiffer evaporation 
spectra and the enhancements of the small ER survival 
probabilities. 

The theoretical understanding of the rapid increase in 
k with A is not clear. Shlomo and Natowitz 65[ assumed 
the effects of long-range correlations wash out when T be- 
comes similar in magnitude to the collective energy Huji of 
each of the modes. For many collective modes, u>i varies 
approximately inversely with the linear dimension, i.e. 
uii ~ A -1 / 3 . Values of k extracted from the predictions 
of Shlomo and Natowitz |65| , shown by the dashed curve 
in Fig. [9] have only a gentle mass dependence and do not 
reproduce our experimental points. 



V. LIGHTER NUCLEI AND YRAST ENERGIES 

Due to the exponential-like dependence of k on mass, it 
seems that the kinetic-energy spectra should described by 
an excitation- independent level-density parameter a e // 
for the lighter nuclei. However, light nuclei have their 
own complications as the spin dependence of E yrast can 
be quite strong. This can cause quite pronounced effects 
on the predicted spectra of a particles which can remove 
appreciable angular momentum from the decaying sys- 
tem. Such effects can in principle be isolated if both 
proton or neutron spectra are also measured as nucleons 
tend to remove very little angular momentum and thus 
are much less sensitive to E yrast . However for the lightest 
nuclei, there is a lot more data available for a particles 
than protons. 
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FIG. 15. (Color online) Inclusive proton and a-particle ki- 
netic spectra in the reaction center-of-mass frame measured at 
angles that highlight compound-nucleus emission. The data 
are associated with 117 Te compound nucleus formed in (a,b) 
14 N + 103 Rh and (c) 40 Ar+ 77 Se reactions. The curves are 
again GEMINlH — \- predictions. The solid curves are the de- 
fault calculations with a distribution of Coulomb barriers, the 
excitation-dependent level-density parameter a e //, and the 
prescription for E yra st{J)- For the short-dashed curves a sin- 
gle Coulomb barrier is used and for the long-dashed curves, 
Sierk's values of E yars t(J) are employed. 



Let us concentrate on the spectra for A— 117 to 59 
compound nuclei in Figs. []~5] to [H GEMINI++ cal- 
culations including the distribution of Coulomb bar- 
rier, Sierk's values of E yarst , and the excitation-energy- 
dependent level-density parameter are indicated by the 
long-dashed curves. Calculations with a constant a e // 
A/7.3 MeV" 1 would be essentially identical to these. For 
protons with minimal angular-momentum effects, one ob- 
tained good agreement with experimental data for the 
117 Te, 106 Cd, and 96 Ru compound systems in Figs. [15] 
to [T7J For the 67 Ga system, the proton spectra are not 
very well reproduced in Figs. [TBT b) and [i"8T c). Actually 
is difficult to understand the evolution of the slope of 
the exponential tails of these proton spectra with excita- 
tion energy within the statistical model. Possibility there 
are experimental problems here or there is contamination 
from other processes. In fact for all these lighter nuclei 
the possibility of contamination exists as the data are all 
inclusive. 
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FIG. 16. (Color online) As for Fig.[l5]but for 106 Cd compound 
nuclei formed in 32 S+ 74 Ge reactions. 
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= FIG. 17. (Color online) As for Fig. [I5]but for 96 Ru compound 



nuclei formed in 32 S+ 64 Ni reactions. 



Consider the 59 Cu data from the 32 S+ 27 A1 reaction 
if Fig. 1191 The evaporation-residue cross section repre- 
sents about 85% of the total reaction cross section at 
E beam =W0 MeV (E*=58 MeV) but decreases to 46% 
at E beam =2U MeV (£'*=11Q MeV) [H]. The remain- 
ing component of the reaction cross section is associated 
with binary-reaction dynamics with various degrees of 
damping and these binary-reaction products evaporate 
protons and a particles [87H89j . Very damped binary and 
fusion-fission reactions are associated with extensive an- 



14 




10 20 30 
e [MeV] 




20 
[MeV] 



FIG. 18. (Color online) As for Fig.llSlbut for 67 Ga compound 
nuclei formed in 40 Ar+ 27 Al and 55 Mn+ 12 C reactions. 



gular distributions and thus light-particle emission from 
these processes will not have a strong angular distribu- 
tion and would be difficult to separate from those associ- 
ated with evaporation residues. Clearly not all the inclu- 
sive a and p spectra can be associated with evaporation 
as is assumed in most analyses. The exact extent of this 
contamination from binary reactions has not been estab- 
lished, but in this work, it will be assumed that it is not 
large for a particles and the basic features of the spectra 
can be traced to evaporation from the fused system. 

For a particles, the GEMINIH — h predictions signifi- 
cantly over estimate the yield in the high-energy tail for 
many of the data sets. In fact these predicted spectra do 
not have exponential tails in the sense that the spectral 
tails decrease linearly on a log plot. This is an indica- 
tion that the predicted enhancement of the high-energy 
region is not a consequence of high temperatures, but of 
angular-momentum effects associated the steep increase 
of Eyrast with J. The angular-momentum effects are 
most pronounced for the more symmetric reactions such 
as the 40 Ar+ 27 Al reactions in Figs. lTSl' a) andfTSJc) which 
populate a region of E* — J space near the yrast line at 
high spins. One also finds the same for the higher-energy 
32 S+ 27 A1 reactions in Fig. HH] 

A large number of previous studies have noted that 
calculations with Sierk's or the RLDM values of E yarst 
are incapable of reproducing a-particle spectra from light 
systems with large angular momentum [l^, [t7| - |20| . lool — 



FIG. 19. (Color online) As for Fig. OH but now for 59 Cu com- 
pound nuclei. For the _E*=60-MeV data, the square and cir- 
cular data points represent the results measured at #; a (,=25° 
and 45°, respectively, while for the 110- MeV data they corre- 
spond to #i a 6=15° and 30°. 



Huizenga et al. [52j reproduced experimental a- 
particle spectra by using a modified yrast energy given 
by 



E. 



yrast 



(J) 



-{l + 5iJ 2 + S 2 J 4 ) 



(21) 



which contains two free parameters, 82 and adjusted 
for each compound nucleus. Equally good fits to the 
data can be obtained by using the Sierk calculations out 
to an angular momentum J* and subsequently allowing 
Eyrast (J) to increase linearly for higher spins, i.e., 



E 



yrast 



(J) 



Esierk(J) if J < J* 

K Esierk(J*) + ( J — J*)E' Sierk {J*) if J > J* 

(22) 

This has the advantage of having only one free parame- 
ter making interpolation and extrapolation easier. Also 
with increasingly large values of </*, the effect turns off 
as Sierk's calculations become more linear (see later) . In 
addition if J* is made larger than the input compound- 
nucleus spin distribution it has not effect. Thus if J» 
increases with A, it allows a smooth transition to heavier 
nucleus where Sierk's values can reproduce experimental 
data. 
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FIG. 20. (Color online) Values of J„, the angular momentum 
for which the Sierk yrast energy is modified, are plotted again 
the mass of the first a-daughter nucleus. The line shows a 
fits to these values which is used in subsequent GEMINlH — h 
calculations. 



The value of J* was obtained from fits to the data from 
59 Cu, 67 Ga, 96 Ru, and 117 Te compound nuclei and the 
values are plotted against the A of the a-daughter system 
in Fig. These data points can be fit with the linear 
function 



J* = 0.319A 



(23) 



shown by the solid line. GEMINlH — h predictions with 
this global parametrization of J* are shown by the solid 
curves in Figs. [TBI to 1191 and reproduce the experimental 
distributions reasonably well. The exception is for the 
106 Cd compound nucleus where the original long-dashed 
calculations in Fig. [IJJ] obtained with Sierk's E yrast (J) 
values produced a better fit. 

In Fig. EU we compare the modified E yra st energies 
to Sierk's calculations for 63 Cu and 55 Co, the daughter 
nuclei following a evaporation from the 67 Ga and 59 Cu 
compound nuclei. In addition are shown values obtained 
by Huizenga et al. obtained from fitting these data with 
Eq. ([2"Tj) [52J. Although the values from this work are 
slightly lower than those of Huizenga et al. at the high 
spins, the most important comparison is that the slopes 
of E yrast (J) are very similar at these high spins. As 
mentioned before, evaporation spectra are not sensitive 
to absolute level density. In this case the calculations 
are sensitive to the J dependence of p which is dictated 
by the spin dependence of E yrast . The E yrast values of 
Huizenga would also give good reproduction of the exper- 
imental data if they were used in GEMINlH — h Evapo- 
ration spectra thus give information on the J dependence 

of E yras i . 

Huizenga et al also suggested that E yrast at these 
high spins not be interpreted as just the rotational-plus- 
deformation energy of the nucleus after shell and pairing 
effects have vanished. Rather they should be treated as 
effective values that may take account of other effects 
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FIG. 21. (Color online) Rotational-plus-deformation energies 
verses the nuclear angular momentum for (a) 63 Cu and (b) 
55 Co. Curves are shown for the dependence caluculated by 
Sierk form his macroscopioc model. These can be compared 
to results obtain from fitting a-particle evaporation spectra 
in this work and by Huizenga et al. 



such as a spin dependence of the level-density parame- 
ter, or spin dependence of collective enhancements. 

In Figs, [lj)] to [TTJ1 the short-dashed curves again show 
the predictions with no barrier distributions in the trans- 
mission coefficients. As the absolute barriers are smaller 
for these lighter nuclei, the effect of the barrier distribu- 
tions on the spectra are reduced. However, the inclusion 
of the distributions (solid curves) still improves the agree- 
ment with the a-particle data except for the 106 Cd com- 
pound nucleus where the long-dashed curves give better 
fits. The 106 Cd data has thus proven exceptional in the 
ingredients necessary to fit both the exponential tails and 
the Coulomb barrier region. The standard calculations 
represented by the solid curves in Fig. [TJJ] would fit much 
better if the experimental spectra were shifted down in 
energy. 

For the 59 Cu and 67 Ga systems, it is clear that the en- 
hancement from the barrier distributions at the largest 
excitation energies and Jo values is not sufficient to re- 
produce the experimental a-particle spectra in Fig. 1191 
Again there are questions about contamination from 
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events not associated with evaporation residues. Majka 
et al. [9(| have investigated the need for a J dependence 
of the transmission coefficients which they associate with 
the increasing deformation of the equilibrium configura- 
tion with spin. At present we have not attempted to 
modify the transmission coefficients as a function of J in 
GEMINI++ to better reproduce the data. 

For the proton spectra, the inclusion of the barrier dis- 
tribution practically has no effect. However, the subbar- 
rier region in the m Te (Fig. H3J), 96 Ru (Fig. [T7J) are still 
underestimated in the calculations even though the rest 
of the spectral shape is well described. Again, there are 
questions as to whether this is a problem with contamina- 
tion from other processes. Alternatively, enhancements 
to the proton subbarrier region can also arise if the evap- 
oration residue is sufficiently proton rich. If the decay 
chain of particles leads to a daughter nucleus with exci- 
tation energy below the neutron separation energy, but 
above the proton value, then subbarrier proton emission 
competes with 7 emission. Such protons are the source 
of the lowest-energy protons in the GEMINI-) — h predic- 
tions for these systems. 

VI. CONCLUSION 

A systematic review of the ingredient necessary to de- 
scribe the shape of proton and a-particle, and some neu- 
tron evaporation spectra was made. In order to describe 
the low-energy yields of the charged particles, transmis- 
sion coefficients associated with a distribution of barri- 
ers were necessary. This was incorporated in a simple 
way into the statistical model assuming a distribution of 
barriers which was assumed to arise from large thermal 



fluctuations. This could include fluctuations in shape, 
density, or surface diffuseness. 

The nuclear level density was described in terms of the 
Fermi-gas formula which is valid for single-particle exci- 
tations. However an effective level-density parameter is 
used that can also account for collective contributions. 
For light nucleus (A <120), the shell-smoothed values 
of a~A/7.3 MeV -1 , obtained from neutron-resonance 
counting at low excitation energies, was also found con- 
sistent with the evaporation spectra. However for heavier 
nuclei at large excitation energies, smaller level-density 
parameters are needed. Evaporation spectra were fit with 
an excitation-energy-dependent level-density parameter 
where the excitation-energy dependence increases very 
rapidly with A. This excitation-energy dependence 
was also found important in understanding the survival 
against fission in very fissile nuclei and allowed reproduc- 
tion of data that previously was thought to require fission 
transients. 

The angular-momentum dependence of the level- 
density is largely defined by the spin-dependence of the 
macroscopic yrast energy. For light compound nuclei at 
large J, modifications to Sierk's and the rotating-liquid- 
drop model values of the rotation-plus-dcformation en- 
ergies which reduces the angular-momentum dependence 
of the level density were needed to describe experimental 
a-particle evaporation spectra. 

These ingredients were incorporated in the GEM- 
INIH — h code to allow a good description of the spectral 
shape of evaporation spectra over all of the periodic ta- 
ble. 

This work was supported by the U.S. Department of 
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